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ABSTRACT 



The energy range of IBIS is a promising ground for 
testing mutual relations of distinct components ex- 
pected in the spectra of high-energy radiation from 



rotation powered pulsars. According to some polar- 
cap models two such components - due to curvature 
and synchrotron emission - may contribute compa- 
rable amounts of power between 15 keV and 10 MeV 
(Rudak & Dyks, 1999). Zhang & Harding (2000) 
argued recently for the inclusion of a third possi- 
ble component, due to inverse Compton scattering 
(ICS) of soft thermal photons on secondary e -pairs. 
Here we present the results of Monte Carlo calcula- 
tions of all three spectral components within a polar- 
cap model which allows for interactions of relativis- 
tic particles with the soft photons coming from the 
pulsar surface. For teragauss pulsars with the sur- 
face temperature of a fewx 10 5 K the ICS component 
dominates the spectrum in the energy range below 
10 MeV, and thus its presence increases the ratio 
of X-ray to 7-ray luminosity (in comparison to the 
models ignoring the ICS on secondary e^-pairs) to a 
level observed in the Vela pulsar. 

Key words: radiation processes: non-thermal; stars: 
neutron; pulsars: general. 



1. INTRODUCTION 



The standard polar cap model of pulsars refers to the 
pair-photon cascades induced by photons of curva- 
ture radiation (CR) due to ultrarelativistic electrons 
leaving the neutron star (NS) surface (primary par- 
ticles). Although successful in reprod ucing gamma- 
ray spectra ob served above 100 MeV ( Daugherty fe 
Harding, 1982) the model fails to reproduce the rel- 
ative level of X-ray and gamma-ray emission. As 
shown by Zhang & Harding (2000) this might be 
caused by neglecting the resonant inverse Compton 
scattering (ICS) of thermal photons from the surface 
by the secondary e pairs. The authors developed an 
approximate analytical method to model the X-ray 



luminosities (with the ICS contribution included) of 
pulsars detected by ROSAT and ASCA. Their results 
indicate that in terms of power in the X-ray range 
the ICS component may be competitive to both the 
curvature and synchrotron compon ents. Howeve r, 



in order t o verify the conclusion of Zhang & Hard- 
ing (200"o| ), and to obtain detailed spectra the Monte 



Carlo simulations are necessary. 

In this work we present basic features of such Monte 
Carlo spectra for different pulsar parameters. In par- 
ticular, the role and properties of the spectral com- 
ponent due to resonant inverse Compton scattering 
on secondary e -pairs are addressed. We find that 
the shape and the strength of this component are 
sensitive to the strength of the local magnetic field 
as well as to the size and temperature of the thermal 
source. 



2. THE MODEL 



As a starting point we took a model of high-energy 
processes in pulsar magneto spheres as proposed by 
Daugherty fc Harding (1982 ). The model was then 



furnished wit h a modern version of a stron g electric 
field E|| after [Harding fc Muslimov (1998[ ) with the 
lower, and the upper boundary of the accelerator at 
the height ho, and h c above the stellar surface, re- 
spectively. In each case the size of acceleration zone 
h c — ho was determined selfconsistently in the course 
of calculations. An extended source of soft thermal 
photons located on the pulsar surface has been added 
(see below for details). Primary electrons slide along 
open dipolar lines of the magnetic field and attain 
ultrarelativistic energies in the presence of this field. 
At the same time they loose the energy via two pro- 
cesses: curvature radiation (CR) and inverse Comp- 
ton scattering (ICS). Due to a large strength of the 
, the former process dominates the cooling of the 
electrons. In consequence, vast majority of electro- 
magnetic cascades is induced by curvature photons 
and very few by ICS (upscattered soft) photons. This 
allows us to ignore any ICS-induced cascades and to 
model the cascades in the following way: 
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First, the CR-induced cascade of e ± pair s and syn- 



c hrotron ph otons is simulated following Rudak & 



Dyks (1999). The primary electrons are accelerated 
along the last open lines of purely dipolar magnetic 
field and emit curvature photons. Some of these pho- 
tons undergo one-photon magnetic absorption pro- 
ducing pairs. The pairs emit in turn synchrotron 
photons (SR) which also may produce higher gener- 
ations of e pairs with their own SR photons, etc.. 

In the second step, the simulation of the resonant 
I CS fo r the pairs proceeds after Dyks & Rudak 
(^000|). Three different cases for the source of ther 



mal photons are considered: A) a hot cap of temper- 
ature Thot — a few x 10 6 K and radius i?hot — 10 5 
cm centered on the magnetic pole; B) entire neutron 
star surface with T sur f ~ a few x 10 5 ; C) the hot cap 
of A) with the rest of the surface at T = 5 • 10 5 K. 
In order to simplify the calculations of scatterings 
in the cases A) and C) (with the hot cap) we made 
the following axial geometry approximation: for the 
local ambient soft photon field as seen by a e^-pair 
we took a field that would be seen if the pair were 
moving along the dipole axis at the same radial co- 
ordinate. 

For the sake of completeness the ICS component due 
to the primary electrons is calculated in the last step, 
even though this component is not important ener- 
getically. 

We would like to emphasize that the numerical code 
constructed for this model does not follow the direc- 
tional pattern of the outgoing radiation. It calcu- 
lates a 'single particle' spectrum of a pulsar, ie. the 
spectrum due to a single primary electron acceler- 
ated to an ultrarelativistic energy in the presence of 
a strong electric field. In this sense our work should 
be considered as a first approximation in modelling 
phase-averaged spectra. 



3. RESULTS 



Fig.l shows the spectrum calculated for a pulsar with 
B pc = 6 TG, P = 0.1 s (the parameters close to 
those of the Vela or B1706— 44) and a hot cap with 
R hot = 10 5 cm and T hot = 5 ■ 10 6 K. The total spec- 
trum (thick solid line) is a sum of three components: 
the curvature spectrum (dot-dashed line), the syn- 
chrotron spectrum emitted by the pairs (short 
HpstipH line), am] the inverse Ckunpion spectrum of 



photons scattered by the e ± pairs (thin solid lino). 



The fourth spectral component - the inverse Comp- 
ton spectrum of photons upscattered by the primary 
electron - shown as the long dashed line is never en- 
ergeticaly important. This is because the primary 
undergoes a strong acceleration and its energy in- 
creases quickly above the value required to fulfill the 
resonant condition. For this reason we have also ig- 
nored any magnetic absorption of these ICS photons 
and the subsequent creation of e^-pairs. At the same 
time the ICS component due to the e pairs is com- 
parable in power to the CR and SR components. The 
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Figure 1. The radiation energy spectrum per loga- 
rithmic energy bandwidth due to a single primary 
electron. The pulsar parameters are B pc = 6 TO 12 G, 
and P = 0.1 s and the case A (R^ot = 10° cm j 
Thot = 5 • 10 6 K) is considered. The spectrum con- 
sists of four components: the curvature radiation 
(dot-dashed) and the ICS (long dashed) are due to 
the primary electron, while the synchrotron radiation 
(short dashed) and the ICS (thin solid) are due to the 
secondary e pairs. Thick solid line indicates super- 
position of the CR, the SR and the e ICS compo- 
nents. The absorbed part of the curvature spectrum 
is marked with dotted line. The dotted verticals mark 
the energy range of IBIS. 



secondary e^- pairs do not accelerate and remain in 
resonance until they lose most of their energy via 
resonant scatterings. 

The e ± -ICS component (thin solid line in Fig. |l|) 
has a well-kno wn shape worked out analytically by 
Dermer (1990] ) for the 'thick target' regime of the 
ICS ic. for the case when scatterings considerably 
decrease the energy of the scattering particle. It re- 
sembles a broken power law with the break located 
at e 1 ^ 3 ~ Bj 2 /T 6 MeV, where B 12 = B/(10 12 G) 
and T 6 = T/(10 6 K). Below the break the pho- 
ton index a ICS assumes the value close to 0, while 
above the break a ICS — —2 (it is the same as for 
the energy distribution of t he IC S emittin g pair s, 
cf. |Daugherty k Harding , |l989| ; [Dcrmer| , |l99C| ). 
The shape of the synchrotr on component is some - 
what more complicated (see Rudak fc Dyks (1999] ) 
for detailed analysis), nevertheless its level lowers 
considerably below the blueshifted cyclotron energy 
get ^ £b!\\ ~ 3 MeV(P/(0.1 s))^Bp C /TG) ( |Rudak 
Dyks| , |1999| ; |Harding fc Daughcrtyl, [19990- Wit h 



the CR break at ~ 10 z MeV ( [Rudak & Dyks, 1999Q , 
most of the power contained in all these spectral com- 
ponents in the case shown in Fig. [I] is confined to the 
energy range above ~ 10 MeV. 

This paucity of power in hard X-ray/soft gamma-ray 
range, clearly visible in Fig.l, stays in disagreement 
with observed spectra of high-i? (ie. classical) pul- 
sars (note the extreme case of the Crab pulsar with 
a comparable energy out put within soft X-ray a nd 
hard gamma-ray window, Thompson et al., 1997). 
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Figure 2. The radiation energy spectra per logarith- 
mic energy bandwidth for three cases of the thermal 
photon source mentioned in the text: A - the hot cap 
with Rhot — 10 5 cm and Thot = 5 • 10 6 K; B - the 
entire surface at T sur f = 5 • 10 5 K; C - the hot cap of 
A with the rest of the surface at T = 5 • 10 5 K. In all 
cases B pc — 6-10 12 G and P = 0.1 s. The dot-dashed 
line is for the curvature component, the dashed line is 
the synchrotron component and the thin solid line de- 
notes the inverse Compton spectrum of photons scat- 
tered by the e^ pairs. Dotted line presents contribu- 
tion of thermal radiation from the surface. Thick 
solid line is the total spectrum per single primary 
electron. 



In order to increase the spectrum level in this range 
we will proceed with three different modifications of 
the case presented in Fig.l. 

First, we will change the size and temperature of the 
soft-photon source, ie. we will consider the cases B 
and C. In Fig. 2 the spectrum of Fig. 1 is presented 
with more details (note the wider energy range) in 
the upper panel and compared with two spectra cal- 
culated for case B (the middle panel) and case C (the 
lower panel). The dotted lines in Fig. 2 present the 
blackbody spectrum of the thermal photon source. 
Its normalization was estimated in an approximate 
way described elsewhere (Rudak & Dyks (2001), in 
preparation). The other lines have the same mean- 
ing as in Fig. 1. Note that the ICS component for 
the warm surface case (thin solid line in the mid- 
dle panel) extends well into the X-ray regime and at 
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Figure 3. The radiation energy spectra per loga- 
rithmic energy bandwidth for the same cases of the 
thermal photon source as in Fig. but for B pc 
{).{') • 10 12 G and P = 0.1 s. The meaning of all lines 
is the same as in Fig. ^. 

10 keV its level exceeds by two orders of magnitude 
a level of the spectrum for the hot cap case (upper 
panel). It has no longer a single power-law shape. 
The significant difference between the shape of the 
ICS component for the hot cap case (upper panel 
in Fig. 2) and the warm surface case (middle panel) 
comes from a difference between the energy loss rates 
of electrons in these two cases. The energy-loss-scale- 
length A-y is at the resonance of the order of 10 2 cm 
and 10 4 cm in the two cases, respectively. In the 
hot cap case (A) electrons lose most energy over a 
very short distance from the creation point and then 
leave the resonance regime. In the warm surface case 
the pairs lose energy in a much slower rate and fulfill 
the resonant condition over the substantially longer 
length scale. Because of the larger size of the ther- 
mal photon source (the entire surface) they can still 
scatter resonantly even at a few stellar radii above 
the surface where the soft tail of the ICS component 
is produced. When both the hot cap and the rest of 
surface emit thermal X-ray photons, then both the 
break at S^ r s and the soft tail are present in the ICS 
component of the spectrum (the case shown in the 
lower panel of Fig. 0). 

Second, we can force the spectrum to extend towards 
lower energies by either lowering the strength of B pc 
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ing level for B pc = 6 • 10 G. We emphasize that the 
main reason for this effect is not a low value of the 
energy loss rate due to the ICS in weak magnetic 
fields (in both cases the pairs lose almost 100% of en- 
ergy due to scatterings) but the smaller total energy 
content of the 7ymc 2 distribution for cascade pairs 
in the low-B case (cf. Zhang fc Harding, 2000). This 



is also the main reason for which the e^-ICS com- 
ponent is negligible for millisecond pulsars (ie. the 
photon index predicted for the INTEGRAL r ange is 
that of synchr otron component: a SR = —1.5, Rudak 



fc Dyks (1999D ) 



Fig. H is to illustrate the case of elevating the accel- 
erator to some high altitude. Three subcases, with 
/io = 1, 2, and 3 R ns were considered. Moreover, 
the figure enables the comparison of the model pre- 
dictions with the data for the Vela pulsar, co ming 
from the compilation by Thompson et al. (1997) and 
Strickman et al. (199S )' The spectra were calcu- 
lated for the entire surface at T sur f = 1.26 • 10 6 K 



( Ogelman, 1995 ). In spite of large altitudes at which 



the pairs are created (low B) , the e ICS compo- 
nent contributes considerably within the X-ray range 
and improves agreement between the theory and the 
data. 
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Figure 4- The radiation energy spectra of the Vela 
pulsar (B pc = 6 • 10 12 G, P = 89 ms, T surf = 1.26 • 
10 6 K). Three panels correspond to different altitudes 
of a lower boundary of acceleration zone: ho = 1, 2, 
and 3 R ns above the surface. The meaning of all 
lines is the same as in Fig. [|. Filled dots (some 
with error bars) present d ata from RXTE, OSSE , 
COMPTEL and EGRE T ([Thompson et al. (199% , 
Strickman et al. (199!\ )). Triangles denote OSSE 



and COMPTEL upper limits. 



at the polar cap or elevating the accelerator to some 
high altitude (ho > 0) to lower the local magnetic 
field strength. The latter case is justified by the idea 
of the high a ltitude accelerator as introd uced and 
argued for by Harding fc Muslimov (1998] ). It is the 



e -ICS component which reacts in this desirable way 
when the strength of B decreases: this component 
extends towards lower photon energies faster than 
the SR component (e]£. s oc B 2 /T whereas e c t oc B). 

Let us lower the magnetic field strength at the polar 
cap by one order of magnitude: B pc = 0.6 ■ 10 12 G. 
The spectra for the cases A, B and C are now 
presented in Fig. ||. In accord with the relation 
f{^ s oc B 2 /T, the ICS component is shifted towards 
the lower photon energy by two orders of magnitude 
and dominates for £ ^ 1 MeV. Again, the break at 
£ br S i s present for the hot cap case (T = 5 • 10 6 K, up- 
per panel) and there is a soft tail and no sign of the 
break for the warm surface case (T = 5 • 10 5 K). The 
level of the e ± -ICS component for B pc = 0.6 • 10 12 G 
is one order of magnitude lower than the correspond- 
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